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Abstract 

We present a complete analysis of K ^ 37r7 decays to 0{p'^) in the low-energy expansion 
of the Standard Model. We employ the notion of "generalized bremsstrahlung" to take full 
advantage of experimental information on the corresponding non-radiative K ^ Sir decays. 
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1 Introduction 



The present experimental status of K ^ 8717 decays is rather meager. So far, only the two 
channels with a charged kaon in the initial state have been detected experimentally with 
very low statistics |^. None of the decay modes of a neutral kaon have been seen. 

This unsatisfactory experimental situation will change soon, especially after the comple- 
tion of the e~^e~ collider DA$NE in Frascati. In this ^-factory one expects Q a total yield 
of 7.5 ■ 10^ KlKs pairs and 1.1 ■ 10^° K^K~ pairs per year. For up-to-date information on 
the future prospects of kaon physics we refer to Ref . [|| . 

With a sufficient number of events, what can one learn from a study of those decays? The 
appropriate framework for such an investigation is chiral perturbation theory (CHPT). 
To lowest order in an expansion in momenta and meson masses, the radiative decays are 
completely determined by the non-radiative amplitudes for — > Stt. At next-to-leading 
order, a full-fledged CHPT calculation of nonleptonic weak amplitudes of 0{p'^) is required 
(cf., e.g., Ref. 0). Among other ingredients to be discussed in Sec. ||, important components 
are the one-loo|) amplitudes with a single vertex from the lowest-order nonleptonic weak 
Lagrangian and tree-level amplitudes due to the corresponding Lagrangian L^^^^~^ 

of O(p^). 

There are three main issues we want to address: 

i. Bremsstrahlung completely determines the lowest-order amplitude, but it also con- 
tributes at next-to-leading order (and at higher orders as well). Is there a unique pro- 
cedure to use all the available information on the non-radiative amplitudes to O(p^), 
either from experiment or from theory? The answer is positive as shown previously for 
a general radiative four-meson process 0. Here, we put the concept of "generalized 
bremsstrahlung" 0] to a practical test. 

ii. The nonleptonic weak Lagrangian of 0{p^) contains a number of low-energy constants 
0, |l^ that are little known at present. Can we expect to extract relevant information 
on those constants from K — >■ 8717 data? 

iii. More generally, can one make definite predictions for these radiative kaon decays within 
the Standard Model? 

The outline of the paper is as follows. In Sec. ^, we set up the kinematics and discuss 
the low-energy expansion of K — » Svr and K 37r7 amplitudes up to 0{j)^). We discuss the 
concept of generalized bremsstrahlung that takes full advantage of the available experimental 
information on the non-radiative amplitude in the form of a fourth-order polynomial in the 
momenta. In Sec. |^, we calculate the electric tree-level amplitude of 0{p'^) in terms of the 
appropriate low-energy constants. We give a fairly complete list of experimentally accessible 
radiative kaon decays that depend on those weak constants of 0{p'^). The calculation of the 
electric loop amplitude is deferred to an Appendix. To the same order in the chiral expansion, 
the magnetic amplitude is a pure tree-level amplitude that receives both direct (local) and 
reducible (nonlocal) contributions. These are put together in Sec. ^. Numerical results for 
rates and spectra of the four transitions occurring at 0{p'^) are collected in Sec. |[ Some 
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conclusions are presented in Sec. p. All relevant formulas for the one-loop amplitudes are 
contained in an Appendix, recapitulating and applying the results of Ref . . 



2 Low— energy expansion 

^ '^i{Pi)'^2{P2)'^3{P3)l{k) is specified by five scalar 
^ = P4{pi-P2), ti = k-pi (2 = 1,..., 4) (2.1) 



The kinematics of the decay K{—p4 
variables which we choose as 



s = {pi+ P2) 



with 



i=l 



i=l 



Any three of the tj together with s and u form a set of independent variables. 

The transition amplitude can be decomposed into an electric and a magnetic part: 

A{K ^ 37r7) = ee^{k){E^ + e^^^M^n (2.2) 

with 



k>'E„ = 



e^,,^k^M''f"^ = . 



To lowest order in the chiral expansion, the amplitudes for both radiative and non- 
radiative transitions are generated at tree level by the effective chiral Lagrangian of 0{p^), 



C2 + C^ 

The strong part has the well-known form H 



A5|=l 



£2 = — {D^UD^U^ + 2BM{U + U^)) 



(2.3) 



(2.4) 



where ( ) denotes the trace in three-dimensional flavour space. F is the pion decay constant 
in the chiral limit (F ^ = 92.4 MeV), M. is the quark mass matrix and B is related to 
the quark condensate. The unitary 3x3 matrix field U incorporates the eight pseudoscalar 
meson fields. In the exponential parametrization. 



U = exp (zV2$/F) 



V2 V6 



$ = $ 



t 



TT 



V 



TT 



K- 



V2 Vq 



K+ \ 



2ri 

76 / 



(2.5) 



with Kl = = {K^ + K^)/^ and Ks = = i{K^ - K^)/\/2 in the limit of CP 
conservation. For the processes under consideration, the covariant derivative D^U can be 
restricted to 

D^U = d^U + ieA^[Q,U] 
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with the photon field and the quark charge matrix Q. 

The weak \AS\ = 1 Lagrangian in (p.3|) can be written in the form (our notation and 
conventions are those of Ref. [§) 

C[^'^=' = GsF'iXL^L^) + 6-27^' (^^^23^^^^! + ^^/.2i^^3) + h.c. , (2.6) 
A = i(A6-a7) , L, = iU^D,U . 



The couphng constants Gg, G27 in (|2.6| ) measure the strength of the octet and the 27-plet 
part, respectively, of the strangeness changing weak interactions. From K titi decays one 
finds 

IGsl ~ 9 ■ 10"^ GeV-2 , ^27/^8 ~ 1/18 . (2.7) 



At lowest order, the magnetic amplitude iVf in ( p^.2D vanishes since there is no e tensor 
in the Lagrangian (|2.3|) . The electric amplitude, on the other hand, is completely determined 
by the corresponding non-radiative amplitude A{s, v) via Low's theorem |TT| : 

+ 2 A12 + (Am - A24) 

+ 0{k) (2.8) 

with (the meson charges in units of e are denoted g^, with I]f=i Qi = 0) 

4 M 

A% = A% = {q,t^-q,U)D^ 



A';- = -^^'. = 7^-^- (2.9) 

Since there are no terms of 0{k) at lowest order in the chiral expansion, the leading-order 
electric amplitude is completely determined by the explicit terms in ( |2.8| ) usually called 
"internal bremsstrahlung" . 

At next-to-leading order, 0{p^), the situation is much more complicated. A nonleptonic 
weak amplitude of 0{p'^) receives in general four types of contributions P]: 

i. Tree-level amplitudes from the effective chiral Lagrangian C[^^^-^ of 0(/) with the 
proper octet and 27-plet transformation properties. 

ii. One-loop amplitudes from diagrams with a single vertex from £2^'^'"^ in the loop. 

iii. Reducible tree-level amplitudes with a single vertex from C^^^^~^ and with a single 
vertex either from the strong Lagrangian £4 or from the anomalous Wess-Zumino- 



Witten Lagrangian |12 
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iv. Reducible one-loop amplitudes, consisting of a strong loop diagram connected to a 
vertex of C^^^^~^ by a single meson line. A typical diagram of this type contains an 
external K — n or K — rj transition, possibly with an additional photon (generalized 
"pole diagrams"). The calculation of such diagrams is simplified by a rediagonalization 
of the kinetic and mass terms of C2 + £2^'^'"^ ("weak rotation" |[T3[|). 



For the decays K —>■ Sttj, all four mechanisms are relevant. Most of them also appear 
in the non-radiative amplitudes. Via Low's theorem (|2.8| ), the non-radiative amplitude of 
C(p^) will contribute to the electric part of the radiative amplitude. Unlike at lowest order, 
this is however not the whole story at 0{p^). The question then is how to use in an optimal 
way the amplitude A{s, v) of O(p^), either from theory or from experiment, for calculating 
the radiative electric amplitude of the same order. 

In a recent paper 0, we have presented the general theoretical framework for the treat- 
ment of radiative four-meson amplitudes like K Sn'j. The essential point is the concept 
of "generalized bremsstrahlung" , 

E>^ = E^^ + 0{k) , (2.10) 
where Eq-q is defined in terms of the non-radiative amplitude A{s, v) |0: 

= A(.,.)S'^ + 2^^A^, + ^^(A5',-A^J 

+ 2^?^^M5'4-tiA^4] • (2.11) 

Referring to Ref. for a more thorough exposition, we concentrate here on the practical 
advantages of generalized bremsstrahlung. Many of the terms in the above list of four mech- 
anisms appear in both the radiative and the non-radiative amplitudes and are therefore 
automatically included in Eq^. This is in particular true for most of the renormalization 
parts that are trivially carried over from A{s, v) to E^, but also for many of the so-called 
reducible contributions (items iii and iv in the above list). For instance, all the weak low- 
energy constants A^,, [|10| contributing to both K — Svr and K 37r7 are completely taken 
into account by Eq^. Therefore, only the genuine radiative low-energy constants A^i4, . . . , 
A^i7 will show up in E^ — Eq^. 

In the following we use the experimental A' — > Stt amplitudes to derive Eq-q. If we had 
limited ourselves to an analysis at the center of the Dalitz plot of A' — > Svr data |T^, |1^, or 



just to linear slopes |T^, there would have been no need to extend (^^ ) to (|2.11|) . However, 



the quadratic slopes are observed and the K ^ Sti amplitudes are written as polynomials of 



second order in s and u [|^, |T9], ^ to fit the experimental data. The second derivatives in 
( p2.1lD are thus needed to take advantage of all the experimental information available from 
K Sir. The (electric) direct emission term E^^ — Eq^ is then a genuine radiative part of 
the amplitude not related to the non-radiative transition. 
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In the numerical analysis we have used the following parametrization of the K 
amplitudes pH, pT 



A{K+ - 
A{Kl- 
A{Ks- 



7lW + ) 



77 Vr TT 



0\ 



with 



vr vr TT 



TT TT TT 



X = 2u/Ml, 



+ c,(y2 + xV3) - (4 - - xV3) , 

2a,(l + iao + + [6,(1 + 2/3o) + 62(1 + iS^W 

+ 2ce(y2 + XV3) + (4 + - XV3) , 

a„(l + iao ~ i«o^) — + ^/3o)V' 
+ CniV + XVS) - 4(1"' - XV3) , 



y = (s - So)/M2, 



So 



(2.12) 



(2.13) 



The numerical values for a,, etc. (in units of 10 ^) are given by [0, |21 



Cc 
62 



-95.39 ±0.40 
24.47 ±0.34 , 
0.68 ±0.17 , 
-1.63 ±0.34 , 
-3.91 ±0.40 , 



On 
bn 

dn 

do 



84.35 ±0.57 , 
-28.11 ±0.49 
-0.05 ±0.22 , 
1.27 ±0.45 , 
0.21 ±0.51 . 



(2.14) 



For the phases associated with the absorptive parts in ( |2.12[ ) we use the lowest-order CHPT 
predictions ao = 0.13, a'^ = -0.12 and /3o = -5o = 0.047 g]. 

The decomposition ( p.l2| ) is based on isospin symmetry. Moreover, the numerical values 
in ( p.l4| ) have been obtained by a fit |19[ where for simplicity the imaginary parts were set 
to zero. Present data on — 37r are too poor]] (especially in the Kg channel) both to relax 
the assumption of isospin conservation and to be sensitive to the small imaginary parts. As 
a consequence, our numerical predictions for the generalized bremsstrahlung amplitudes in 
K S'K'-f are affected by systematic errors and must be considered as preliminary. A new 
detailed analysis should be performed when complete and accurate K ^ Stt data will be 
available. 

In the next two sections we discuss separately the electric and the magnetic amplitudes 
for the various channels. To 0{p'^), there are four non-vanishing transitions: 



Kl 



TT 



K+ 
Ks 



71 71 71 7 



TT 



We make the following simplifications for the calculation. The 27-plet part of the nonleptonic 
weak Lagrangian is not included in the calculation of direct emission amplitudes, i.e. in 
[E^^ — Eq-q) and M^''" . This is an excellent approximation in view of the A/ = 1/2 rule. 
Moreover, in the loop diagrams we have only kept the dominant two-pion intermediate 
states. Since the loop amplitudes will turn out to be rather small anyway, this restriction is 
justified a posteriori. Finally, CP conservation will be assumed throughout the analysis. 

^ Note that we have not taken into account the very recent and accurate results of Serpukhov-167 |^ in 
the K'^ 7r°7r°7r+ channel. 
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3 Electric amplitudes 



To O(p^), the electric amplitude can be written as 

E'^ = Eq^ + -E^ountor + -^loop,subtracted • (3-1) 

Use of the generalized bremsstrahlung amplitude Eq-^ greatly simplifies the calculation of 
both the tree-level and the loop part of (|3.1| ). For instance, all the reducible contributions 
(items iii and iv in the list of Sec. ^ to the electric amplitude are automatically contained in 
Eq-q. This can be shown almost without any calculation by going back to the definition ( p.ll ) 
of generalized bremsstrahlung. The only exception that needs some (tree-level) calculations 
are amplitudes proportional to the strong low-energy constant Lg [^] with an external weak 
transition. Although there are strong radiative four-meson amplitudes proportional to Lg, 
the explicit calculation shows that they do not contribute to K Stt'-/ after a weak rotation. 

Another consequence of using generalized bremsstrahlung in ( p.l[ ) is a much simpler 
form of £^^ounter- ^11 the low-energy constants appearing in both radiative and non-radiative 
amplitudes are already contained in Eq-^. Therefore, only the genuine radiative terms in the 
octet Lagrangian of C(p'^) fl^ 

£f ^1=1 = GsF^ + h.c. , (3.2) 

i 

with dimensionless coupling constants Ni and octet operators Wi, contribute to ^^counter- 
In particular, going through the Lagrangian (|3.2|) one finds that only the four low-energy 



constants A^i4, . . . , can occur in -E^ounter- The relevant parts of the Lagrangian (|3.2|) are 
listed below. 

With Ff^u = d^Ay — duA^ the electromagnetic field strength tensor, the explicit coupling 
for 7r^7r*^7r+7 is given by 

- '-^{Nu - N,, - - N,j)F^,K+n^d^'n^d''7r- . (3.3) 
The corresponding expression for K"^ 7r+7r+7r^7 reads 



(N,, - N,, - N,, - N,j)F^,K+n-d''n+d''n- . (3.4) 



The decay Kl —>■ n^n^ir^'j receives a contribution from 

- ^(A^M - ATis - N^, - N,j)F^,KLid''n%- d" n+ - 27r°9^7r+a'^7r-) , (3.5) 
and Ks 7r+7r^7r'^7 from 

- ^[7{N,, - N,,) + 5{N,, + N,,)]F^,Ksd^n\n-d^Ti+ + ti+O'^ti-) . (3.6) 



In order to facilitate the comparison with other radiative kaon decays, we list in Table |I| 
the combinations of low-energy constants Ni governing the various experimentally accessible 
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Table 1: Kaon decay modes to which the couphng constants Ni contribute. For the Svr final 
states, only the single photon channels are listed. For the neutral modes, the letters L or S* in 
brackets distinguish between Kl and Ks in the limit of CP conservation. 7* denotes a lepton 
pair in the final state. If a decay mode appears more than once there are different Lorentz 
structures in the amplitude. The combinations with are scale dependent compensating 
the scale dependence of the corresponding loop amplitude. The other combinations are scale 
independent. 



TT 


27r 


3n 




7r+7* 
ttS* (S) 
7r+77 


7r+7r°7* 
ttOttS* (L) 

Tl+TT^'-f'-f 

7]-+7r^77 (5*) 
7]-+7r*^7 

7r+7r"7* (L) 
7r+7r"7* (S) 
7r+7r°7* 


7r^7r+7r^7 
7r+7r'^7r°7 
7r+7r^7r'^7 (L) 
7r+7r^7r°7 (S) 


^[4 - Nl, 

14 ' io 
)7 

A^i4 - A^is - A^ie - A^iy 

J7 
5) 

7(iV[,-iV[e) + 5(iV[5 + iVi7) 
iV[4 - Nl, - 3(iV[g - iViv) 
iV[4 - iV[5 - 3(iV[e + iViy) 
iV[4 + 2iV[5-3(iV[6-iVi7) 




7]-+7r~7 (L) 
7]-+7r'^7 


7r+7r~7r°7 (S) 
7r^7r+7r~7 

7r+7r^7r°7 {S) 
7r+7r^7r'^7 (L) 


iV29 + A^31 

3Ar29 - ATgo 
5Ar29 - ATgo + 2Ar3i 

6iV28 + 3iV29 - 5iV3o 



channels. This Table is a slightly extended version of the one appearing in Ref. |^. As one 
can see from the Table, the specific combination of coupling constants A^i4 — A^i5 — Ni^ — Nn 
occurs also in the amplitudes for — > 7r"'"7r°7 and Ks — > 7r"'"7r~7. On the other hand, 
7(A''[4 — A''[g) + 5(A''[5 + A''i7) is a characteristic combination^ for Ks 7r+7r^7r°7 only. 

Both combinations are not yet known phenomenologically. To get a feeling for the typical 
size of these couplings one may appeal to the factorization model that predicts [|T^ 

[Ari4 - ATis - ATie - ATi^]™ = -kf^ = -7 ■ , (3.7) 

[7(Ari4 - A^ie) + 5(Ari5 + ATi^)]™ = 4%^, (3.8) 

where /cj is a fudge factor which naive factorization sets equal to one. Note the potentially 
large counterterm amplitude in Ks — > 7r+7r^7r'^7. Table |I]also indicates that the combination 
NiA - A^i5 - A^i6 - ^^17 is scale independent while 7(A^[4 - iV[g) + ^{N{^ + Nn) is not. 
Consequently, the loop amplitudes are all finite for K^ 7r'^7r°7r+7, K^ 7r+7r+7r^7 and 

^We remind the reader that Nn is scale independent px| |. 
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Figure 1: One-loop diagram for a four-meson transition. For the radiative amplitude, the 
photon must be appended to every charged meson line and to every vertex with at least two 
charged fields. For the case of K ^ 37r(7), —pa = —p4 is the kaon momentum, the other 
three being the pion momenta. The weak (strong) vertex Vi (V2) is defined in Eq. ( |A.1| ), 
with the appropriate coefficients for the various diagrams given in Table 

Kl 7r+7r^7r°7, but divergent for Ks 7r+7r^7r°7. This divergence is renormalized by the 
counterterm combination 7{Ni4 — Niq) + 5{Ni5 + Nn). In the limit where 6*27 is set to zero, 
the two-pion loop does not contribute to the Ks decay . Since we have not included the other 
loop contributions that are numerically negligible, the amplitude of C(p^) for the Ks decay 
is superficially scale dependent. We shall come back in Sec. ^to investigate numerically the 
effect of this scale dependence. 

Finally, the loop contributions to (|3.1| ) have to be calculated. Once again, many contri- 
butions are already contained in Eq-q. The only type of diagram that has to be calculated 
explicitly is shown in Fig. |I| where a photon can be appended to all (charged) lines and 
vertices. In this diagram, Vi is a weak vertex from C^^^^~^ and V2 is a strong vertex from 
£2- Of course, such diagrams without a photon contribute also to the K ^ 3n amplitudes 
of O(p^). In accordance with the definition of generalized bremsstrahlung in (|2.11| ), the ap- 
propriate part has be subtracted from the radiative loop amplitude to obtain -Eioop.subtractcd 
in the complete amplitude (|3.1|). 

The calculation of the loop amplitudes is rather involved in the radiative case. We have 
given in Ref. [0 a compact expression for the radiative loop amplitude with general vertices 
Vi, V2 of 0{p^). In an Appendix, we reproduce the main steps for arriving at the final 
amplitude, together with the relevant vertices for K 37r7. 



4 Magnetic amplitudes 

The magnetic amplitude in (|2.2|) receives contributions from direct and reducible diagrams 



[|^, corresponding to type i and iii, respectively, in the classification of Sec. |[ 

The direct parts (type i) are generated by the operators • • • , in (^.21). Their 
contribution to K^ 7r'^7r'^7r^7 is given by 

- ^(3iV29 - N,,)F^,d^K-'d''n-Ti\' , (4.1) 
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Figure 2: Reducible diagram contributing to the magnetic amplitude at 0{p'^). A weak cubic 
vertex of and an anomalous vertex with three mesons and a photon are connected by 

a single meson line. 



where = e^upaF^" (£^0123 = +!)• The corresponding expression for — >• tt+tt+tt 7 
reads 

- ^(iV29 + N,,)F,,d^K+ d'^Ti^Ti-Ti- . (4.2) 
The decay 7r+7r^7r°7 receives a contribution from 



F2 

and i^'s ^ 7r+7r~7r''7 from 
2ieG: 



-(6iV28 + 3iV29 - 5iV3o)F^.9^i^L5"7r%-7r+ , (4.3) 



F2 



F^,Ks[{hN2s, - Nso + 2N;^)d^Tx\- d"" 7r+ - 2(iV29 + N^x^d^i^^ d"" n'] . (4.4) 



Following the theoretical arguments given in the coupling constants in the anomalous 
parity sector of 0{p'^) can be estimated as 



Aran _ ^ Aran _ ^ 



Q (4-5) 

Aran _ Aran _ ^ 



where the dimensionless coefficients are expected to be positive and of order one. 

The second class of diagrams contributing to the magnetic amplitude are the reducible 
ones (type iii). These amplitudes are due to diagrams with a single meson line between a 
weak I AS* I = 1 vertex and an anomalous vertex from the Wess-Zumino-Witten (WZW) 
functional [1^. For K — > Sn'j, all such diagrams have the structure shown in Fig. ^: a weak 
cubic vertex and an anomalous vertex with three mesons and a photon. 

In the case of ti^ti^ti^'j, there is only one reducible contribution at 0{p'^): the 

kaon emits a neutral and a charged pion, where the tt~^ subsequently makes an anomalous 
transition to 7r''7r+7, 

j-^-L- weak 0/ + WZW + \ 

> TT [7T^ > TT 7r^7) . 

The corresponding amplitude is local because the vr'^vr^ vertex vanishes on-shell 

(remember that we are setting G27 = in direct emission amplitudes). Thus, the complete 
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magnetic amplitude (adding the direct term generated by ([4.1|) ) takes the form 

M^P'^{K+ ^ 7r°7r°7r+7) = ^^(3^2 - - 2)k''plpl . (4.6) 

There are two types of reducible diagrams contributing to ti^tt^tt^'j: the can 

make a weak transition into a real tt^ and a virtual 7r° (or 77) which is then transformed into 
a TT+vr^ pair and a photon, 

weak -I-/ n WZW -I- — \ 

> 7r^(7r > TT^TT 7) , 

T^A- weak _L / WZW ^ — \ 

A ^ TT (?7 y Tl TX 7) . 

The total magnetic amplitude is now given by 

M^''-(ir+ ^ 7r+7r+7r-7) = ,^k^Ma, + 2a,)pl 

HM', - M^)(^I— + -^^)] , (4.7) 



with 

si4 = {Pi + Vif = iy + h + {Ml + - s)/2 , 
S24 = iP2 + Pif = -y + h + {Ml + - s)/2 . 

For fC^L ^ 7r^7r^7r°7 one may either contract the anomalous KiKs'Tf^'y vertex with the 
weak KsTi^Ti^ vertex, or the weak Kl —>■ tt^tt^ transition with the 7r^7r^7r°7 WZW vertex: 

WZW / weak -L _ \ 

Kl — >iT -f{Ks — > TT^vr ) , 

weak J- / — WZW —On 

Kl > TV [TT > TT TT 7j , 

weak _ / + WZW + \ 

Kl > TT (tT > TV TT 7j . 

The last two diagrams give again a local amplitude for a similar reason^ as for K~^ — >■ 7r''7r'^7r+7 
in ([4.6|). Together with the contribution from (|4.3|), we arrive at the magnetic amplitude 



M'^'"^{Kl - TT+TT-TT^) = ^^[24ai + 3a2 - SOaa - 2 - ^^^^ ^f- ^ ^pK . (4.8) 



Finally, we turn to ^ tt+tt n^'-f. In this case, the reducible diagrams have the following 
structure: 

weak + / — WZW — \ 

A5 > TV [TT > TT TT 7) , 

weak — / + WZW + \ 

Ks ^ TT (tT > TT TT 7) , 

weak n / n WZW + _ \ 

Ks > TT [71 > TT^TT 7) , 

weak / WZW + — \ 

As' > TT (77 s> TT TT 7) . 

•^The on-shell amplitude for Kl — » tt+tt^ vanishes in the limit of CP conservation. 
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Combined with (|4.4|) , we obtain 

M'^^-(ir5 TT+TT- A) = ^^A;'^{(5a2-6a3 + 4a4-2)(p2-pi)^P^ (4.9) 

+4(Mi - M,^)(^i— - 

where 

S34 = {PZ + Pif = S + 2{ti+t2) . 



5 Numerical results 

Our numerical results for the various channels are displayed in Tables 0-0- The first column 
shows the photon energy range. In the second column, the contribution to the decay width 
generated by the generalized bremsstrahlung amplitude £'gb in is listed, together with 



the corresponding errors due to the uncertainties of the if — Svr parameters in (|2.14|) . The 
next column shows the relative change of the result if only the Low amplitude (|2.8| ) is used 
instead of Eq-q- In the fourth column we see the effect of adding the electric counterterms 
(using kf = 1 m. ( p.7| ) and ( |3.8| )) and the residual pion-loop contributions -Eioop,subtracted in 
( |3.1|) . Fm in the next column denotes the contribution to the decay width from the magnetic 
amplitudes (for = 1); there is no interference between electric and magnetic amplitudes 
as long as the phase space integration is performed "symmetrically" . 

For the branching ratios in the last column we distinguish between the three channels 
where the leading-order amplitude is not suppressed and the decay Ks 7r"'"7r~7r°7 with a 
suppressed bremsstrahlung amplitude. In the first group of transitions, the dominant 0{E.y) 
effect is given by the difference Fgb — Flow. This deviation from Low's theorem, i.e. from a 
pure QED prediction, could possibly be observed in the near future. In the above channels, 
the residual pion-loop contribution suffers from relatively large theoretical uncertainties: the 
smallness of phase space amplifies isospin-breaking effects generated by the mass difference 
MttO — M^+. However, the effect of -Eioop.subtracted is always so small that it can hardly be 
detected. The contribution of -Ecounter, evaluated within the factorization model, is of the 
same order as -Eioop.subtracted. For — > 7r°7r°7r+7 there is an almost complete destructive 
interference between loops and electric counterterms, while for 7r+7r+7r^7 we find 

^'counter ^ ^loop.subtracted- Finally, in the Kl channel ^ioop,subtracted is bigger than ^counter for 
large E-y. For small E^, the two amplitudes are comparable. Probably only large deviations 
from the naive expectation kf ~ 0{1) could be observed. Also the magnetic contribution is 
very much suppressed in these channels: the ratio Fm/Fgb is typically smaller than 10~^. 

Interference effects between electric and magnetic amplitudes could in principle be larger. 
For instance, observables like det(pi,p2,P3,P4) (for the decays with three different pions in 
the final state) or u det{pi, p2, P3, Pa) (in the case of two identical particles tti, 112) are sensitive 
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to such interferences. To 0{p^), the interference term 

^.^p<x(^^M,V + E;^V) (5.1) 

is proportional to the relatively small absorptive part of the electric amplitude. Thus, the 
leading-order piece of i?^ does not contribute in ( p.l[ ). Nevertheless, the possibility of inter- 
ference measurements should be kept in mind once sufficiently high statistics will have been 
achieved. 

For the three channels under consideration, the amplitude is completely dominated by 
generalized bremsstrahlung. In the last column of Tables |^ - ^ we therefore list the branching 
ratios based on generalized bremsstrahlung only, corresponding to Fgb in the second column. 
The contributions to the branching ratios from direct emission are completely concealed by 
the present experimental uncertainties of the K ^ 3n parameters. 

Within those errors, our predictions are consistent with standard bremsstrahlung and 
with the available experimental results. Our theoretical branching ratio for K~ — >■ 7r°7r°7r~7 
in Table ^ for E^^ > lOMeV can be compared directly with the experimental result |]l| 

BR{K- 7T%\--f) = (7.4+^1) • 10"^ , E^>W MeV . (5.2) 

For — s> 7r+7r+7r~7, Barmin et al. have reported the branching ratio 

BR{K+ ^ 7r+7r+7r-7) = (1.10 ± 0.48) ■ 10~^ , > 5 MeV , (5.3) 

to be compared with our theoretical prediction 

BR(K+ ^ 7r+7r+7r-7)|theor. = (1-26 ± 0.01) • 10"^ , E^ > 5 MeV , (5.4) 

whereas Stamer et al. ^ have found 

BR(K+ ^ 7r+7r+7r"7) = (1.0 ± 0.4) ■ 10"^ , E^ > 11 MeV . (5.5) 



For the decay Ks 7r+7r^7r°7 the situation is quite different. To lowest chiral order, 
the amplitude can only proceed through a Al = 3/2 transition (via bremsstrahlung) and 
is therefore suppressed by the Al = 1/2 rule. Consequently, the next-to-leading order 
contributions generated by octet operators are becoming relatively more important^. At the 
one-loop level, two-pion intermediate states do not contribute. Therefore, the 0{p^) part of 
the electric amplitude is essentially determined by the counterterm 

iVx,(/i) := - N'[,) + 5(iV[5 + N^r)] (/i) (5.6) 

that is predicted to be large by the factorization model in (|3.8|) . Its rather modest scale 
dependence, 

3 

NksM = NKsif^i) + ^ln(/ii//i2) , (5.7) 
similar phenomenon occurs in the 7r+7r°7 decay |2j, [27| . 
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Table 2: Numerical results for the decay 7r°7r*^7r+7. The photon energy E-^ and the 

decay widths Fgb, Fm are given in MeV. 



Ej 


Fgb 


Fgb — Flow 
F 


Fe — Fgb 
F 


Fm 


BR 


10-20 


(1.38 ±0.02) ■ 10-19 


1.4 ■ 10-3 


2.2 ■ 10-5 


2.8 ■ 10-25 


(2.60 ±0.03) - 10-6 


20-30 


(4.29 ± 0.06) ■ 10-20 


4.5 ■ 10-3 


8.2 ■ 10-5 


7.5 ■ 10-25 


(8.05 ±0.01) - 10-7 


30-40 


(1.45 ±0.03) ■ 10-20 


9.8 ■ 10-3 


2.5 ■ 10-4 


1.1-10-24 


(2.72 ±0.05) - 10-7 


40-50 


(4.48 ± 0.09) ■ 10-21 


1.8 • 10-2 


2.2 • 10-4 


1.1-10-24 


(8.42 ±0.18) - 10-8 


50-60 


(1.09 ±0.03) ■ 10-21 


2.9 ■ 10-2 


-1.0 ■ 10-3 


6.8 - 10-25 


(2.05 ±0.05) - 10-8 


60-70 


(1.49 ±0.05) ■ 10-22 


4.3 • 10-2 


-6.8 • 10-3 


2.0 • 10-25 


(2.81 ±0.09) - 10-9 


70-80 


(3.48 ±0.12) ■ 10-24 


5.6 ■ 10-2 


-1.9 ■ 10-2 


8.9 - 10-27 


(6.55 ±0.23) - 10-11 


10-80 


(2.01 ±0.03) ■ 10-19 


3.3 • 10-3 


4.5 ■ 10-5 


4.1 - 10-24 


(3.78 ±0.05) - 10-6 



Table 3: Numerical results for the decay tt+tt+tt 7. 



E-y 


Fgb 


Fgb — Flow 
F 


Fe — Fqb 
F 


Fm 


BR 


10-20 


(2.32 ±0.02) - 10-18 


-1.7-10-3 


-4.2 - 10-4 


1.3 - 10-24 


(4.36 ± 0.04) - 10-5 


20-30 


(7.63 ±0.07) - 10-19 


-4.8 - 10-3 


-1.2 - 10-3 


3.2 - 10-24 


(1.43 ±0.01) - 10-5 


30-40 


(2.62 ±0.03) - 10-19 


-9.2 - 10-3 


-2.4 - 10-3 


4.1 - 10-24 


(4.93 ± 0.05) - 10-6 


40-50 


(7.66 ± 0.08) - 10-20 


-1.5 - 10-2 


-4.1 - 10-3 


3.2 - 10-24 


(1.44 ±0.01) - 10-6 


50-60 


(1.43 ±0.02) - 10-20 


-2.1 - 10-2 


-6.2 - 10-3 


1.3 - 10-24 


(2.69 ± 0.03) - 10-7 


60-70 


(7.23 ±0.09) - 10-22 


-2.8 - 10-2 


-8.5 - 10-3 


1.2 - 10-25 


(1.36 ±0.02) - 10-8 


10-70 


(3.44 ±0.03) - 10-18 


-3.4 - 10-3 


-8.5 - 10-4 


1.3 - 10-23 


(6.46 ±0.06) - 10-5 



is compensated by loop graphs with kaon intermediate states which we have neglected. The 
uncertainty induced by this scale dependence, 

iV,,,(0.5GeV)-iV.,(lGeV)^^^^^^ ^^^^ 

^Ks 

with A^™ given by ( |3.8| ) for kf = 1, is certainly smaller than the intrinsic uncertainty of the 
factorization hypothesis. 

The corresponding numerical results are displayed in Table |[ The numbers for Fe are 
obtained from the sum of Eqb (using the central values of the input parameters 62 and d2 
in (|2.14| )) and the aforementioned 0{p'^) counterterm amplitude. Note that the interference 



is destructive and especially pronounced at large values of E^. The contribution of the 
magnetic amplitude is again shown for = 1. For this channel we list the total branching 
ratio BR = (Fe ± Tuj/^totiKs) for the various photon energy bins. We do not give errors 
for these branching ratios because, unlike for the other three channels, the direct emission 
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Table 4: Numerical results for the decay Kl 



TT + TT 7r°7. 



E 


-L GB 


1'gb — 1'low 

r 


1'e — 1'gb 
F 


i M 




10-20 


(1.32 ±0.02) ■ 10-^s 


-4.5 ■ 10-3 


-4.2 ■ 10-4 


2.1 ■ 10-26 


(1.04 ±0.02) ■ 10-4 


20-30 


(4.89 ± 0.07) ■ 10-^9 


-1.3 ■ 10-2 


-1.1 ■ 10-3 


4.4 ■ 10-26 


(3.84 ±0.06) ■ 10-5 


30-40 


(1.98 ±0.03) ■ 10-^9 


-2.5 ■ 10-2 


-1.8 ■ 10-3 


5.0 ■ 10-26 


(1.55 ±0.03) ■ 10-5 


40-50 


(7.33 ±0.11) ■ 10-20 


-4.0 ■ 10-2 


-1.7-10-3 


3.7- 10-26 


(5.76 ±0.10) ■ 10-6 


50-60 


(2.13 ±0.04) ■ 10-20 


-5.8 ■ 10-2 


2.5 • 10-3 


1.7- 10-26 


(1.67 ±0.03) ■ 10-6 


60-70 


(3.39 ±0.06) ■ 10-21 


-7.9 ■ 10-2 


2.7- 10-2 


3.5-10-27 


(2.67 ±0.05) ■ 10-7 


70-80 


(8.04 ±0.15) ■ 10-23 


-9.9 ■ 10-2 


2.2 ■ 10-^ 


9.5 ■ 10-29 


(6.32 ±0.13) ■ 10-9 


10-80 


(2.11 ±0.03) ■ 10-1^ 


-1.0 • 10-2 


-6.6 • 10-4 


1.7-10-25 


(1.65 ±0.03) • 10-4 



Table 5: Numerical results for the decay Ks 



TT TT 





Tgb 


1'gb — 1'low 

r 






BR 


10-20 


(1.29 ±0.34) ■ 10-21 


1.2 ■ 10-2 


1.1-10-21 


6.5 - 10-25 


1.5 - 10-10 


20-30 


(5.15 ±1.28) • 10-22 


3.8-10-2 


3.4 - 10-22 


1.6 - 10-24 


4.7-10-11 


30-40 


(2.34 ±0.53) ■ 10-22 


7.7-10-2 


9.7-10-23 


2.3 - 10-24 


1.3 - 10-11 


40-50 


(9.97 ±2.12) ■ 10-23 


1.2 - 10-1 


2.0 - 10-23 


2.1 - 10-24 


2.9 - 10-12 


50-60 


(3.34 ±0.68) ■ 10-23 


1.6 - 10-1 


2.2 - 10-24 


1.2 - 10-24 


4.6 - 10-13 


60-70 


(6.09 ± 1.22) ■ 10-24 


2.1 - 10-1 


2.3 - 10-25 


3.4 - 10-25 


7.8 - 10-14 


70-80 


(1.62 ±0.32) ■ 10-25 


2.4- 10-1 


1.7-10-26 


1.2 - 10-26 


4.0 - 10-15 


10-80 


(2.18 ±0.55) ■ 10-21 


3.3-10-2 


1.6 - 10-21 


8.2 - 10-24 


2.2 - 10-10 



amplitude matters with unknown theoretical uncertainties (factorization model). 

Remembering the projected DA$NE yield of 7.5 x IQ'^KlKs pairs per year, the Ks —>■ 
tc'^tt'tc^'j decay rate is still too small for the coming generation of kaon experiments. With 
an additional improvement of statistics, some information might be achieved via time- 
interference measurements p8| [K^^s ~^ 7i^7i^7T^'~f) similar to those recently performed in 
the non-radiative case [^, |SD|. Then interference effects between electric and magnetic am- 
plitudes could in principle be measured since a term like 



e;m,p.) 



(5.9) 



is generated. In contrast to j ^A. 



this term is proportional to the leading-order piece of 
Efj^. We stress that even fixed-target experiments, through regeneration, can perform time- 
interference measurements and in this case a larger statistics is expected. Thus, the Ks —>■ 
7i^7i~n^'j decay mode may still turn out to be a valuable probe for kaon physics parameters 
that is not drowned by bremsstrahlung. 
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6 Conclusions 



Anticipating substantial improvements in the statistics of K ^ 8717 decays in the near fu- 
ture, we have performed a comprehensive and complete analysis of these decays to 0{p'^) 
in the low-energy expansion of the Standard Model. To lowest order, 0{p'^), the decay am- 
plitudes are determined by the corresponding non-radiative amplitudes via Low's theorem 
(bremsstrahlung) . At next-to-leading order, there are different contributions to both elec- 
tric and magnetic parts of the amplitudes: loops and tree-level (counterterm) amplitudes, 
reducible and irreducible contributions. 

A major aspect of our analysis is the concept of "generalized bremsstrahlung" that trans- 
fers the available theoretical or experimental information on X — > Stt decays to the corre- 
sponding radiative amplitudes in an optimal way at the level of 0{p^). For the numerical 
analysis, we have used the factorization hypothesis to estimate the relevant low-energy con- 
stants. 

Returning to the three issues addressed in the introduction, we may summarize our 
findings as follows: 

i. In all three channels where the leading-order amplitudes are not suppressed {K^ — > 
7r°7r°7r+7, K'^ 7r+7r+7r~7, — * tt+tt'tt^), generalized bremsstrahlung completely 
dominates the amplitudes to 0{p'^). The differences to the QED prediction (standard 
or internal bremsstrahlung) could be experimentally observed in the forthcoming round 
of kaon experiments, at least from the statistical point of view. 

ii. For the same channels, it will hardly be possible to extract the appropriate combi- 
nations of low-energy constants from experiment in the near future. This conclusion 
hinges, of course, on the assumption that the factorization estimates are not off by 
an order-of-magnitude in amplitude. In contrast, the counterterm amplitude is impor- 
tant for Ks 7r+7r^7r°7, especially if the rather large factorization estimate is reliable. 
However, for this decay mode the branching ratio is probably too small to be detected 
soon. 

iii. As a general conclusion, the Standard Model allows for quite definite predictions for 
radiative kaon decays into three pions. Especially for 7r°7r*^7r+7, K'^ tt^tt'^tt"'-/ 
and Kl 7r"'"7r~7r°7, the accuracy of these predictions is at the moment only limited 
by the precision with which the parameters of the non-radiative decay amplitudes 
are known. For Kg — > 7r"'"7r~7r'^7, there is some theoretical uncertainty related to the 
relevant low-energy constants. 

As soon as more accurate data will lead to better precision for the K ^ Sir parameters, 
the predictions of the radiative amplitudes can be improved accordingly. Although we have 
only considered total rates and photon energy spectra in this analysis, the investigation of 
more subtle effects like the interference between electric and magnetic amplitudes may then 
become feasible. 
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Appendix: Loop amplitudes 

In this Appendix, we collect the main results of Ref. for the calculation of loop amplitudes 
corresponding to the diagram in Fig. |I|. 

First, we calculate the loop amplitude for the non-radiative process K 3tt. In our 
case, —pa = —p4, is the kaon momentum and Vi, V2 are nonleptonic weak and strong vertices, 
respectively. The pion momenta are generically denoted Pb,Pc,Pd- 

We characterize the vertices Vi, V2 in momentum space by constants Oj, 6j: 

Vi = ao + aiPa-Pb + a2Pa-x + asix'^ - M^) + ttiiy'^ - My) + a:i{pl- M^) + ae^pl - M^) 
V2 = bo + hp,-pd + b2Pc-x + h{x^ - Ml) + 64(2/' - M^) + h^{pl - Ml) + hivl - Ml) . 

(A.l) 

With P = pc + Pd, the non-radiative loop amplitude of Fig. |l| can be represented in the 
following form (all external lines are on-shell): 

F{P) = A{M^)[aihpa-Pb + a4biPc-pd + a4h{P^ + Ml - M^) + aoh + a4bo] 
+ A{My)[aibsPa-Pb + a2bsPa-P + asbip^-p^ + a3b2Pc-P 

+aMP^ - Ml + M^) + aobs + a^bo] 
+ 5(P^ M^, My)[aobo + aobiPc-Pd + aiboPa-Pb + aibiPa-PbPc'Pd] 
+ 5i(P^ M^, My)[aob2Pc-P + a2boPa-P + aib2Pa-PbPc- P + aibiPc-PdPaP] 
+ a2b2\pa-PcB2^{P\M^,My)+pa-Ppc-PB22{P\M^,My)] . (A.2) 

The various functions in ([A.2|) are as defined conventionally (in d dimensions): 

I r d^X 1 

A{M) -- 



{B,B,P„g,^B2o + P,P.B22) = \ J ^^^^^ _ ^^^^^^^"^ _ . (A.3) 

We have chosen to express F{P) in terms of the scalar products 

Pa-Pb, Pc-Pd, P"^, Pa-P, Pc-P, Pa-Pc (A.4) 

instead of using kinematical relations to express all scalar products in terms of the two inde- 
pendent scalar variables s, v. Note that the analytically non-trivial part of ( |A.2| ), involving 



the various B functions, contains only the on-shell couplings oq, cti, 02, &0; ^ij The off-shell 
couplings 03,04,63,64 appear only together with the divergent constants A{M). Since these 
terms are polynomials in the momenta of at most degree two, they will enter in the radiative 
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Table 6: Coefficients of the vertices Vi,V2 defined in ( |A.1| ) for the various loop diagrams. 
Only the relevant on-shell coefficients are listed. 



K(-Pa) - 


— S> 7T(ph) + 














'k(x\'k(v\ — 




On 


CLi 










K+ - 


















— i> TT^TT^ 




-2 


-2 


2M2 

TT 


2 


-2 


K+ - 


^ 7r+ + 


















-2Ml 


-2 


-2 
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2 





K+ - 
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^ TT- + 
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-2 
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-2 
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. ttO + 














n^TT^ - 




MI/V2 
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2Ml 

TT 
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-2 



amplitude only through internal bremsstrahlung and will therefore eventually be absorbed 
in -Eqb- The on-shell coefficients for the various channels are listed in Table ^. 

We now turn to the radiative loop amplitude and decompose it into two parts: 



-^loop 



(A.5) 



The amplitude can be expressed through derivatives of the non-radiative loop amplitude 
F in ( A. 2 ) with respect to the various scalar products ( A. 4 ). In some of the following terms, 
the momentum P has to be replaced by P + /c, leaving all scalar products unchanged that 
do not contain P explicitly: 



F{P)W + 
dF 

+ 



F{P + k) - F{P) 



dF 



d{Pa-Pb) 



ab 



d{pa-P) 



d{pc-Pd) 



d{pc-P) 



+ 



qatc 



dF dF 

{P)-qcta^ AP + k) 



I, 



diPa-Pc] 
tat 



d{Pa-Pc 



ac 



d^F 



d{pa-P)d{p,-P) 



iP)D'^. 



d^F 



aP 



d{pa-P)d{p,-P) 



[P + k)Di:p 



(A.6) 
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We have used the definitions ( |2.9D . When P appears as an index (e.g., in K'^p or -D^p), the 
corresponding momentum and charge in ( p.9|) are P and Qc + Qdi respectively. 

The second part of the loop amplitude (|A.5|) cannot be expressed in terms of F or 
derivatives thereof. For the relevant case of equal loop masses (M^. = My = M^^ =: M), 
takes on the following compact form: 

H^' = a2{tbp':-tapt){iq.-qy){2bo + 2bip,-pd + b2Pc-P)C7oiP^-k-P) 

+ h2{q. + qy)[-2p,-PC;^{P\ -k-P) + 2^,G^2(P^ -k-P) -p^PC^oiP', -k-P)]} 
+ h2{tdP''c - tcP'"d){{qx - qy)[2ao + 2a^Pa-Pb + a2{Pa-P + ta)]C^Q{{P + kf, k-P) 
+ aafe + qy)[-2{p,-P + QC^^^^P + k)\ k-P) - 2t,ar2((P + A;)^, A;-P) 
- (Pa-P + ta)C^o((P + fc-P)]} . (A.7) 



The functions Cjj are defined as 



Cij{u,v) - Cij{u,0) 



Q^{u, v) = (A.8) 
in terms of the three-propagator one-loop functions Cij{p^, k-p) for /c^ = 0: 



i J (27r)^ (x2 - M2) [{x + p)2 - M2] [(x + A;)2 - M^] 
= {C2o{p^,k-p)g^,y + C3l(p^ k-p){p^gyp + p^gf^p + Ppfi-^^) 

+ C'32(p^ k-p){k^g^p + k^g^p + fep^f^,,) + ...}. (A.9) 

We recall the following observations from Ref. 0: 

i. The amplitudes in ( |A.(j| ) and in ( [A.?] ) are separately gauge invariant. 



ii. The amplitude is finite and at least of 0{k). It only contains the on-shell couplings 
ao, ai, 02, &o, ^1, i>2 defined in ( |A.1|) and the charges g^, of the particles in the loop. 



iii. The generalized bremsstrahlung part of the loop amplitude is contained in G^. De- 
noting by i?QB (loop) the result obtained by inserting for A{s,v) the on-shell loop 
amplitude (|A.2| ) in Eq. (|2.11| ), the difference 

= G'^-E^bOoop) (A. 10) 

is at least of 0{k). Moreover, by construction of Eq^ all the divergences in A'^ are 
renormalized by counterterms with an explicit field strength tensor. Finally, A^ is 
finite for 02^2 = 0. 

Putting everything together, the subtracted loop amphtude E^^^^ subtracted ( P.IQ is given 

by 

Kop,subtracted = E + ' (A-H) 

loops 

The sum extends over the various configurations listed in Table ^ 
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